if( !require("BiocManager", quietly = TRUE))
  install.packages("BiocManager") #判断包是否安装，没安装先安装包
if( !require("org.Hs.eg.db", quietly = TRUE))
  BiocManager::install("org.Hs.eg.db")
if( !require("clusterProfiler", quietly = TRUE))
  BiocManager::install("clusterProfiler")
load("D:/杨熊冰/R语言/实验报告/class5_volcano.RData")
library(org.Hs.eg.db)
library(clusterProfiler)
prostat$ID[1:5] #查看ID信息,  通过网上查询可知ID为SYMBOL信息
ID <- prostat$ID[prostat$P <0.05] #提取待富集ID
ego <- enrichGO(ID, "org.Hs.eg.db", 
                qvalueCutoff = 1, #由于本数据富集结果不显著,为了后续画图修改阈值
                pAdjustMethod = "none", #由于本数据富集结果不显著,为了后续画图修改校正方法
                keyType = "SYMBOL", #设置ID类型为SYMBOL
                ont = "BP", #选择注释类型为BP
                universe = prostat$ID #设置背景基因
)
dotplot(ego)
ggsave(filename = "富集BP图.jpg")